Perioral secretions enable complex social signaling in African mole-rats (genus Fukomys)

Subterranean common mole-rats of the genus Fukomys (family Bathyergidae) live in large, cooperatively-breeding families. Odor cues have been hypothesized to play an important role in mediating social behaviors in the underground ecotope, but only little is known about the role of olfactory signaling in burrowing mammals. Here we characterize the so far neglected perioral glands of Fukomys and other African mole-rats as an important source of olfactory social information. Histology demonstrates these structures to be derived sebaceous glands that are developed regardless of sex and reproductive status. However, gland activity is higher in Fukomys males, leading to sexually dimorphic patterns of stain and clotting of the facial pelage. Behavioral assays revealed that conspecifics prefer male but not female perioral swabs over scent samples from the back fur and that male sebum causes similar attraction as anogenital scent, a known source of social information in Fukomys. Finally, we assessed volatile compounds in the perioral sebum of the giant mole-rat (Fukomys mechowii) via GCxGC-MS-based metabolomic profiling. Volatiles display pronounced sex-specific signatures but also allow to differentiate between intrasexual reproductive status groups. These different lines of evidence suggest that mole-rat perioral glands provide complex odor signals which play a crucial role in social communication.

www.nature.com/scientificreports/ While there is no evidence that Fukomys use perioral secretions to actively scent mark their environment, they could nevertheless play an important role in social communication. Both Cryptomys and Fukomys engage in conspicuous reciprocal cheek nuzzling in various social situations (Fig. 1E). This behavior often initiates copulation (see e.g., 28 -Cryptomys hottentotus; 29 -Fukomys mechowii) but is also observed when unfamiliar individuals, regardless of sex, meet for the first time ( 30 -Fukomys anselli) and when family members greet each other after short periods of separation (pers. obs.). In the mating context, it typically precedes anogenital inspection 29 . Obviously, perioral scents are of interest to both partners during these interactions, suggesting a role in social and particularly sexual signaling.
Several studies hypothesized that only certain status groups of mole-rats display visible perioral secretions. In particular, it has been proposed that just reproductive males and thus one animal per family, would exhibit perioral stains. Although this idea has never been tested explicitly, various field studies have relied on this character to identify breeding males (F. anselli- 31 ; F. damarensis-32 ; F. mechowii- 33 ). The restriction of this feature to breeding males would suggest a role in both intra-and intersexual signaling and could imply an involvement in sexual suppression of subordinate males. However, there is no consensus about whether perioral secretions are indeed specific to male breeders. For instance, Kawalika 34 reported that in Zambian giant mole-rats (F. mechowii) perioral stains are displayed by both sexes irrespective of reproductive status. On the other hand, Caspar et al. 11 noted anecdotally that the degree of expression in perioral secretions of captive Ansell's mole-rats (F. anselli) is sex-specific but unrelated to breeding status, with males in general displaying more intense stains than females (see also 35 for F. mechowii). In any case, an exclusiveness of secretions to particular status groups in mole-rat communities would have important implications for their social function.
Here, we aim to characterize the occurrence, composition, and biological significance of perioral secretions in mole-rats of the genus Fukomys by the aid of morphological and histological observations, behavioral assays, and chemical analyses.

Materials and methods
All statistics were performed in R 36 .

Histology of the mouth corner integument in African mole-rats (Fukomys spp. and
Heterocephalus glaber). We sampled the perioral integument of the mouth corners in 13 Fukomys individuals, comprising three species (F. anselli, F. mechowii, and F. micklemi) as well as both sexes and intrasexual status groups (breeder vs. non-breeder), to gain insights about the morphology of glands producing perioral secretions in these animals (see Suppl. Table 1, for further data on sampled specimens). For comparison, we also included samples from four naked mole-rats (Heterocephalus glaber), another cooperatively-breeding species of bathyergid in which perioral glands but no visible secretions have been reported 24,27 . All animals were adults deriving from the Department of General Zoology in Essen and were sacrificed for other research projects (respective animals were decapitated in deep ketamine/xylazine anaesthesia-compare 37 ). No animals were sacrificed primarily to obtain perioral samples. Hence, the representation of species, sexes and reproductive status groups is imbalanced. Ultimately, we recovered perioral glands in the tissue samples of all of the four examined individuals of F. mechowii and H. glaber, respectively, in three out of six F. anselli and in none of the three F. micklemi (Suppl. Table 1). Given that we recovered glands across status groups and sexes in both genera, we are convinced that the apparent lack of glands in some individuals reflects issues with tissue sampling rather than their absence in the respective animals.
Skin was excised and prepared for standard histological sectioning and staining. Fukomys mouth angles were shaved with a handheld trimmer (Isis GT420; Aesculap, Suhl, Germany) before sampling. Tissue samples were fixed in 4% buffered paraformaldehyde for 24 h at 4 °C, subsequently transferred to 1× DPBS buffer (PAN-Biotech; Aidenbach, Germany), and stored at the same temperature until being automatically dehydrated (Tissue Processor TP12-RWW Medizintechnik; Hallerndorf, Germany) and embedded into paraffin (EG1150 H embedder-Leica Biosystems; Deer Park, USA). Embedded samples were manually sectioned (thickness: 5 µm) on a Microm HM 340 rotary microtome (Microm International; Walldorf, Germany), transferred to a digital precise water bath (Witeg WB-11; Wertheim am Main, Germany) warmed to 40 °C and subsequently mounted on glass slides. Tissue sections were stained using a standard hematoxylin-eosin manual staining protocol. ROTI ® Histokitt II (Carl Roth; Karlsruhe, Germany) was used as a xylol-based cover medium for the tissue sections. Samples were examined and photographed on a VHX-600 digital light microscope (Keyence; Osaka, Japan).
We used ImageJ 38 to take quantitative measurements of glandular cell sizes from the micrographs. We measured the area of medially sectioned mature non-pyknotic cells in perioral glands from Ansell's mole-rats (n males = 2, n females = 1) and naked mole-rats (n males = 1, n females = 3). Subsequently, we tested for sex differences as well as species differences in this variable. For Ansell's mole-rats, we also checked for differences in cell size between perioral glands and ordinary, hair follicle-associated sebaceous glands. We compared cell sizes between gland types by means of the two-sided t-test and calculated Cohen's d as a measure of effect size. To explore effects of species and sex on perioral gland cell size we computed a generalized mixed effect model using the lmer() function from the lme4 package ( 39 ; version 1.1.31; gamma distribution) of the following form: cell size ~ species + sex: species. The individual ID was included as a random factor. Normality of data as well as of model residuals was checked with the Shapiro-Wilk test.
Occurrence of perioral staining among sexes and status groups of Fukomys. The  www.nature.com/scientificreports/ (F. micklemi) was studied to test the influence of selected biological variables. The two species represent distantly related congeneric lineages 40 . We examined mole-rats with monitored life histories living in the laboratories of the Department of General Zoology, University of Duisburg-Essen, Essen, and the Department of Zoology, University of South Bohemia, České Budějovice (Suppl. Table 2). All mole-rats were housed in social groups with food provided ad libitum.
The giant mole-rats descend from animals caught in the Zambian Ndola region and exhibit the diagnostic karyotype of 2n = 40. We included 30 males (14 reproductive ones, 16 non-reproductive ones) and 68 females (14 reproductive ones, 54 non-reproductive ones) of giant mole-rats, resulting in a total sample of n = 98. The imbalance among the sexes and the two female status groups is a result of the strongly female-biased sex-ratio of neonates found in this species 35 , which also affected our sampling efforts for the mass spectrometric analyses (see below). The Micklem's mole-rat lab lineage derives from animals caught at Kataba in western Zambia, the type locality of the species 41 , and are characterized by a karyotype of 2n = 60. In this species, we studied 40 males (19 reproductive ones, 21 non-reproductive ones) and 32 females (20 reproductive ones, 12 non-reproductive ones), thus comprising a total sample of n = 72.
For documentation of the perioral stains, animals were briefly separated from their group, weighed, and photographed. Based on these photographs the degree of expression was scored on a species-specific qualitative scale from 1 (no visible secretion) to 4 (excessive secretion). Classifying criteria are enumerated in Table 1 and stain categories are visualized in Fig. 2.
For each of the two species, we separately calculated cumulative link models for ordinal regression (clm() function in the ordinal package-42 ; version 2022.11.16; logit link function) to estimate the effects of biological variables on the expression of perioral stains. We used a two-step approach to the models: First, we used reproductive status, sex and the interactions between these two factors as model predictors to answer the question whether perioral stains are a sex-specific status signal (model I: stain pattern ~ sex * sex : reproductive status). Subsequently, we tested whether body mass (in g, log-transformed) and individual age (in months, logtransformed) predicts this trait within the sexes (model II: stain pattern ~ sex : log 10 (age) + sex: log 10 (body mass)). www.nature.com/scientificreports/ Olfactory preference tests. We ran olfactory preference assays to assess the relative informative value of mole-rat perioral secretions compared to other body odours. Adult (at least 17 months old) Micklem's mole-rats (F. micklemi) participated in this part of the study, because the species expresses notable perioral stains while being of small size and thus easily deployable for testing. The olfactory preference assay was designed as a twochoice set up, in which test subjects were presented with odorous swabs taken from different body regions of the same foreign conspecific donor animal with visible perioral secretions (see below). One option was invariably constituted by swabs taken from the perioral area, while the second one either derived from the donors' dorsal pelage or perineal area. While dorsum samples acted as a simple control to simulate the presence of a foreign conspecific, anogenital smear is known to convey complex social information in Fukomys 13,16,43 .
Odorous smear from the donor animals was collected with moistened commercial cotton swabs that were gently rubbed against the respective body region. For perioral and dorsum samples, this was done until the tip became visibly stained. During the procedure, one experimenter would briefly fixate the donor animals while a second one collected the samples. The swabs were rolled out on the surface of glass cuvette lids (10.2 cm × 8.5 cm × 1.1 cm) which were subsequently used to present the odors to the test subjects (analogous to 13 ). During this procedure as well as during the set-up of the two-choice assay, the respective experimenters wore gloves to avoid olfactory contamination of the equipment.
Test subjects were individually taken from their home cages and brought to a darkroom in which the assay set up was deployed. To allow the experimenter to operate, the room was illuminated by an LED table lamp emitting monochromatic red light (Parathom R50 80.337 E14 Red 617 nm, 6 W, Osram; Munich, Germany), which is invisible to African mole-rats 44 . The animals were placed in a terrarium (50 cm × 38 cm × 30 cm) in which they were presented with the two glass plates carrying the odors of the donor animal. Glass plates were positioned equidistant from the center along the long-axis of the terrarium with a 5 cm distance to the walls and were fixed with tape on the underside to remain in place. The position (left vs. right) at which the different odor types were presented was randomized. Test animals were observed exploring the set-up for three minutes after being placed into the center of the terrarium. All experiments were recorded (SONY HDR-CX505 camcorder) and behaviors were quantified based on these recordings. Interest in the presented odors was approximated by the time spent sniffing at the respective glass plate. Sniffing was defined as the animal lowering and moving its head over the glass plate accompanied by visible movement of the rhinarium. Besides sniffing time, the latency until first sniffing for either glass plate and the number of sniffing events was quantified. Later on, however, these measures were deemed to be uninformative and not considered for further analyses, since the animals were for the most part alternating between the two presented options. Experimental runs in which total sniffing time was < 5 s were discarded, leaving us with 66 valid runs in total (Suppl. Table 3). The researcher quantifying the sniffing responses was blind regarding the identity of the offered scent samples.
Animals were tested in three situations. The sample sizes itemized for sex are shown in brackets: 1. Dorsal vs. perioral secretion, male donor (n males = 13, n females = 15): Mole-rats could choose between swabs from the dorsal pelage and perioral region of a foreign male conspecific. 2. Dorsal vs. perioral secretion, female donor (n males = 12, n females = 12): Mole-rats could choose between swabs from the dorsal pelage and perioral region of a foreign female conspecific. 3. Anogenital vs. perioral secretion, male donor (n males = 13, n females = 11): Mole-rats could choose between swabs from the perineum and perioral region of a foreign male conspecific. Males rather than female donors were chosen because greater interest in male secretions was found in previous runs comparing perioral and dorsal samples.
Deviations in the sample compositions for the different test situations derive from changes in the lab population caused by deaths and animals being transferred to other institutions. If possible, individual animals were tested across all three situations. A maximum of one experimental run per animal per day was performed. After each run, the set up was cleaned with water and mild detergent.
Differences in sniffing time for perioral compared to dorsum and anogenital samples were statistically assessed for each of the three test situations. Data were checked for normality using the Shapiro-Wilk test. Parametric datasets were analyzed using the paired t-test and by calculating Cohen's d as a measure of effect size, nonparametric ones by aid of the paired Wilcoxon signed rank test and Wilcoxon r to indicate effect sizes. Responses of males and females were compared for all test situations. However, sex differences were not found to be significant and thus data from males and females were pooled for all analyses to increase statistical power. Total sniffing times were compared across the three test situations by means of the Kruskal-Wallis test.
Metabolomic profiling. We collected perioral secretions from giant mole-rats (F. mechowii) to identify volatile organic compounds which might facilitate social communication via two-dimensional comprehensive gas chromatography with mass detection (GCxGC-MS). Samples from 26 animals were analyzed (Suppl. Table 4; secretions from four further animals were used to calibrate the GCxGC-MS). F. mechowii was selected for this aspect of the study since secretions are expressed in greater quantities compared to congeneric species. Furthermore, we applied the same methodology to analyze volatiles from small amounts of hay and cereals to consider potential diet-related contamination of secretions.
Samples were collected from living, manually restrained animals. Perioral secretions glue the hair in the cheek region together, so that clotted hairs could be swiftly cut and manually collected in Eppendorf tubes. Subsequently, samples were stored at − 20 °C until analysis. The dynamic headspace method was used to sample the secretion compounds. The sampling process was carried out automatically using a multi-purpose sampler device (MPS, Gerstel, Germany). Sebum samples were placed in 10 ml glass vials and incubated for 5 min at 50 °C www.nature.com/scientificreports/ before a flow of nitrogen of 20 ml/min was used for continuous volatile extraction. The extraction was carried out for 10 min. Volatiles were sorbed on a Tenax sorbent packed in a glass tube (Tenax® TA, Gerstel, Germany) at 20 °C and subsequently released in a thermal desorption unit (TDU) at 295 °C into a programed temperature vaporizer (PTV) inlet precooled to a − 30 °C where the volatiles were captured on a glass wool. The PTV inlet unit was then fast heated up to 300 °C and the analytes were introduced into a gas chromatograph. Volatiles were then analyzed employing the GCxGC-MS (Pegasus 4D, Leco Corporation, USA). A combination of nonpolar and polar separation columns was used for the separation: Primary column: Rxi-5sil MS (28 m × 0.25 mm ID, Restek, Australia); Secondary column BPX-50 (1.6 m × 0.1 mm ID, SGE, Australia). Other parameters were set as follows: splitless mode, constant He flow 1 ml/min, modulation time 3 s (hot pulse 0.9 s), modulation temperature offset with respect to the secondary oven 15 °C. The temperature program applied on the primary oven: 35 °C (hold 1 min), then increase (8 °C/min) to 320 °C (hold 2 min). The temperature offset applied on the secondary column was + 5 °C. Transferline temperature was held at 250 °C. The mass detector was equipped with an EI ion source and TOF analyzer enabling a unit mass resolution. The scanned mass range was 29-700 m/z. The ion source chamber was held at 280 °C. LECO's ChromaTOF v4.5 was employed to control the instrument and for data processing. Selected compounds were identified by automatically matching their mass spectra with a library of mass spectra (NIST MS 2.2, USA).
To prepare the bioinformatic analysis of the GCxGC-MS data, we first generated histograms of all samples and blanks. The resulting distribution was bi-modal with compounds that occurred only in samples (green line in Fig. 6A) and those that occurred in both samples and blanks (intersection). To decide which compounds were true positive metabolites, we used the mixtools package ( 45 ; version 1.2.0) which calculates the posterior probability (p < 0.05) for the identity to either of the two peaks within the mixture of two overlapping normal distributions. Next, we applied a normalization based upon quantiles, which normalizes a matrix of peak areas (i.e. intensities) with the function normalize.quantiles of the preprocessCore package ( 46 ; version 1.60.1). We used this type of normalization because the Gaussian filtering resulted in a clear bell-shaped curve and there were no global differences in the distributions between the groups caused by biological variation. To explore potential sources of variation in our data, we used sparse partial least squares discriminant analysis (sPLS-DA) within the mixOmics package ( 47 ; version 6.22.0) for the fact that it has satisfying predictive performances with large datasets. To extract p-values of differentially abundant compounds, we used the power law global error model (PLGEM-48 ) which is an efficient tool to calculate differentiation within large data sets (e.g., proteomes, transcriptomes, metabolomes) with distributions that deviate from normality (see methods in 49 ) for more details). We used ggplot2 ( 50 ; version 3.4.0) to visualize differentially abundant compounds.
Ethics statement. Ethical review and approval for the behavioral assays was not required since animal housing (approved by permit no. 32-2-1180-71/328 Veterinary Office of the City of Essen) as well as all experiments complied with the corresponding animal testing regulations and were approved by the animal welfare officer in charge. All behavioral tests as well as handling protocols conformed to the relevant ethical standards and did not harm the animals. All applied methods are reported in accordance with the ARRIVE guidelines.

Results
Histology of the cheek region in African mole-rats. We found large, specialized sebaceous glands in the mouth corner regions of both sexes and irrespective of reproductive status in Fukomys (Fig. 3A-E) as well as Heterocephalus (Suppl. Fig. 1). The exocrine glands in question are strongly branched, multilobated structures (Fig. 3A). Their wide secretory ducts open directly onto the skin surface ( Fig. 3B; Suppl. Fig. 1). Although we were unable to measure their full extent, these glands form fields covering an area of at least several square millimeters in sections from both species. As all sebaceous glands, they show a holocrine secretion pattern, releasing lysed cell masses into their glandular ducts (Fig. 3C). When stained using HE-solution, the glands appear well demarcated and in a light violet color. Individual gland cells are large (see below), contain well visible nuclei, and increase in size from their formation site in the peripheral layer to the center of a lobule. Perioral glands are embedded in fibrous connective tissue permeated by striated muscle cells. Otherwise, the histology of the cheek region was unremarkable and our observations aligned with those of earlier studies on bathyergid skin 25,27,51,52 . Regular sebaceous glands (Fig. 3F), but no other types of skin glands, were also recovered in both genera. In Fukomys, most hairs are arranged in follicle compounds, which are associated with one to several small, globular sebaceous glands (compare 52 ). In Heterocephalus, body hair is almost completely reduced so that regular sebaceous glands are only found in association with the vibrissae 25 .
Apart from deviations in morphology, cell size between perioral and regular sebaceous glands differed significantly in F. anselli. Cells from regular sebaceous glands had a median area of 135.4 µm 2 (n = 69; SD: 33.71) in the studied sections, while it was 314.7 µm 2 (n = 279; SD: 92.16) for perioral gland cells (t = 26.105, p < 10 -15 ; d = 2.23). There were no significant differences between the sexes in the size of the cells constituting ordinary sebaceous glands (t = 0.218, p = 0.831; d = 0.09).
Results of the linear mixed effect model on perioral gland cell sizes in bathyergids are summarized in Table 2. The regression model revealed perioral gland cell size to be smaller in Heterocephalus compared to Fukomys (t = − 2.824, p = 0.005). We also found a significant sex difference in cell size in Heterocephalus (t = 4.055, p < 0.001), with the male exhibiting larger cells (n = 92; median: 365.5 µm 2 ; SD: 113.76) than the females (n = 222; median: 192.1 µm 2 ; SD: 72.76). In Fukomys, no such dimorphism was recovered (t = 0.100, p = 0.921).
Occurrence of perioral staining among sexes and status groups of Fukomys. We found perioral stains to be highly sex-specific but not related to reproductive status in both studied Fukomys species. Patterns www.nature.com/scientificreports/ of sex and status-dependent perioral stain expressionsare shown in Fig. 4 and are itemized in Supplementary  Table 3. Males tended to show a greater development of stains than females and excessive perioral secretions (category 4) were exclusively found among the males of both species (Fig. 4). Secretion was found to be exaggerated and more strongly sexually dimorphic in F. mechowii compared to F. micklemi (Fig. 4).
The results of the ordinal regression analyses on stain expression are provided in Table 3. In our initial models, we found that sex is a good to excellent predictor of perioral stain intensity in both species, with males showing a more pronounced expression than females (Table 3, p < 0.02). However, reproductive status does not influence stains in either sex or species (Table 3, p > 0.2). The subsequent models tested whether individual age or body mass might intrasexually influence the expression of stains. We found that these factors showed no significant effects in F. micklemi, regardless of sex, while we found age to be a significant predictor of stain expression in male F. mechowii (p = 0.02; Table 3).
Note that we only sampled adult individuals here and did not systematically study when stains formed during ontogeny. Our anecdotal observations suggest that stains manifest at an age between 12 and 18 months. The youngest individual in which we noticed perioral stain was a 7-month-old F. mechowii female, which was sampled for the characterization of volatile compounds (Suppl. Table 4).  www.nature.com/scientificreports/  www.nature.com/scientificreports/ Olfactory preference tests. Mole-rats showed great interest in male but not female perioral swabs, although individual differences in responses were pronounced (Fig. 5). The animals showed a highly significant preference (paired Wilcoxon test; V = 353, p < 0.0001, r = 0.76) for odor derived from male perioral secretions (median sniffing time: 23.17 s; SD: 26.3) over swabs from the dorsum (median sniffing time: 12.45 s; SD: 14.9), irrespective of their sex. There was one dropout run for this condition (final sample: n males = 13, n females = 14). We recovered a pronounced difference in the median sniffing time of female test animals (43.6 s) compared to males (14.1 s) for the perioral swabs, but it nevertheless failed to be significant (Wilcoxon test; W = 63, p = 0.19, r = 0.26). In fact, significant sex differences were recovered in none of the three experimental conditions. In contrast to male-derived samples, interest in female perioral swabs (median sniffing time: 10.25 s; SD: 13.0) was not different from that in swabs taken from the dorsal pelage (median sniffing time: 11.79 s; SD: 14.6) of donor animals (Wilcoxon test; V = 58, p = 0.39, r = 0.21). Additionally, we noticed the largest number of dropouts for this condition, with 6 out of the total 10 dropouts being observed here. This left us with 18 valid runs (n males = 7, n females = 11). Given this sex-specific difference, we continued with testing preferences for male perioral secretions compared to anogenital smear, which is a known source of social information in mole-rats. There was no significant difference in sniffing times between male perioral swabs (median sniffing time: 24.1 s; SD: 15.24) and the perineal swabs (median sniffing time: 18.05 s; SD: 15.19) examined by the mole-rats (t-Test; t = 1.38, p = 0.18, d = 0.24). Three dropout runs occurred in this condition (final sample: n males = 10, n females = 11).
The Kruskal-Wallis test did not indicate significant differences in the total time spent sniffing at the two offered odor samples, across the three test situations (Kruskal-Wallis Χ 2 = 3.93, p = 0.14).
Metabolomic profiling. GCxGC-MS profiling of perioral secretions in F. mechowii yielded a total of 765 volatile compounds. However, some metabolites were also partially present in blanks. So, we employed Gaussian  Table 5).
To explore whether GCxGC-MS profiles reflect sex, we employed a supervised form of discriminant analysis, sPLS-DA, that relies on the class membership of each observation. The sPLS-DA model, based on significant components, accounted for 15% (component 1) and 9% (component 2) of the data variance (Fig. 6B). We used the Area Under the Curve (AUC) to provide evidence that the discrimination is excellent in both dimensions (AUC1 vs. AUC2) thus yielding AUC = 0.916 and p = 0.0004 for component 1, and AUC = 1, p = 0.00002 for component 2. This analysis thus indicates a strong sexual dimorphism in perioral sebum volatiles. Our AUC approach also revealed that component 2 is more informative than component 1 when looking at separations based on the reproductive status within each sex (Fig. 6C). Female non-breeders are most different from the remaining status groups (AUC = 0.99, p = 0.0001) but all others could also be reliably differentiated (female breeders-AUC = 0.96, p = 0.004; male breeders-AUC = 0.88, p = 0.0025; male non-breeders-AUC = 0.88, p = 0.037). It should be pointed out that male non-breeders cluster completely within the odor space of breeders, while the differentiation among female status groups appears far more pronounced (Fig. 6C). Yet, each separation is significant on component 2 and thus volatile profiles may have the potential to signal reproductive status in both sexes.
To further test the hypothesis that males and females have different profiles of volatiles we used PLGEM models of differential expression/abundance to extract levels of sexual dimorphism. A volcano plot (Fig. 6D) visualizes the striking differences between the sexes, which are already detectable at the bottom of this highly symmetrical plot. However, strictly statistically speaking, when fold difference is set to 2 and α = 0.05, a total www.nature.com/scientificreports/ of just 28 compounds are sexually dimorphic with 7 compounds being male-biased and 21 ones being femalebiased. Next, we asked whether those sexually dimorphic compounds are also highly abundant. Thus, we recalculated signal intensities to abundances (0-100%). In Fig. 6E, we clearly see that those compounds that are most abundant (> 0.29%) are least likely to be sex-biased, or in other words, sexual dimorphism is expressed by many compounds with rather low abundances while the species-level odor-space comprises non-dimorphic compounds which are highly abundant. The comparison of volatiles between perioral sebum, hay, and cereals demonstrated a similar number of identified compounds in the food items and the secretions (Suppl. Fig. 2; Suppl. Table 6). Yet only ~ 10% (n = 68) of the total compound diversity was shared between all the three sets and even fewer being exclusively shared between sebum and hay (n = 2) as well as sebum and cereals (n = 3). We can thus exclude notable biases due to food intake for our chemical analysis.

Discussion
Histology of perioral glands. This study is the first to report on the presence and histology of enlarged and complex perioral glands in male as well as female Fukomys. For Heterocephalus, they had been previously discussed in passing by Quay 24 and Kimani 27 . The anatomy of the perioral glands complies to a pattern reported from all hystricomorph rodent taxa studied so far (Capromys pilorides, Hystrix indica, Myocastor coypus 24,25 ), in that exclusively sebaceous and no sudoriferous gland components constitute the structures. However, to our best knowledge, secretions of perioral glands in all these taxa, or in fact any other rodent, do not permanently clot the fur in a way similar to Fukomys, suggesting that condition to be exceptional. Interestingly, the mouth corners appear to be devoid of comparable secretions in Heterocephalus.
The finding of sexually dimorphic perioral gland cell sizes in Heterocephalus was surprising, as this species displays remarkable monomorphism in other morphological and physiological traits (e.g., 53 ). Sex differences in gland cell size could indicate that specific secretion variables might deviate between male and female Heterocephalus. Yet, similar differences are not evident in Fukomys, although gland activity is markedly higher in males than in females of this genus (Fig. 4, see below). In any case, our findings must be interpreted with caution, as the number of sampled individuals is very low for both genera and only includes a single Heterocephalus male as well as one Fukomys female. Although Heterocephalus has been extensively studied in regard to its social organization and communication 54 , chemical signaling in these animals remains essentially unknown, so that perioral secretions might represent a promising subject for future research.
Apart from the conspicuous perioral glands we describe here, no integumental scent glands have been formally documented in Fukomys so far. However, it is known that various hystricomorph rodents, including Heterocephalus, possess specialized sebaceous skin glands in the perianal area 25,27 . Tullberg 55 even described such anal glands in the Southern common mole-rats of the genus Cryptomys (= "Georychus coecutiens"), providing additional indication for their presence in Fukomys. Secretions from these glands might well underlie the great social significance of anogenital scents in Fukomys 13,16 and can be expected to complement olfactory signaling via perioral secretions.
Occurrence of perioral staining among sexes and status groups of Fukomys. We found perioral stains to be strongly sexually dimorphic but not affected by reproductive status in either of the two studied Fukomys species. Thus, although both sexes possess well developed perioral glands, the quantity of secretion is typically far greater in males. This might suggest that sexual hormones affect perioral secretion patterns. Indeed, it has been shown that the activity of sebaceous glands in rodents is stimulated by androgens and inhibited by estrogens 56 . For instance, androgens have been demonstrated to increase the size and activity of the supracaudal sebaceous gland of the guinea pig (Cavia porcellus), another hystricomorph species 57 . The specialized perioral sebaceous glands of Fukomys might similarly respond to these hormones, giving rise to the observed sex differences in fur stain at the mouth corners.
Our analyses demonstrate that reproductive status does not notably influence the expression of perioral stain in adults of either sex. Hence, we challenge the assumption that this is a characteristic trait of breeding males 8,31-33 . But if adult male breeders and non-breeders essentially display equally noticeable perioral secretions, why do field studies often report it from just one animal per colony? A simple explanation might lay in the dispersal behavior of wild Fukomys. At the time when perioral secretions start to be notably developed, which appears to be at an age between 12 and 18 months, many male non-breeders have already dispersed from their natal family. Long-term field studies on F. damarensis indicate that male dispersal happens at a mean age of 12 months already, thus limiting the time that a non-breeder with fully developed perioral secretions might be captured in its natal colony 5 . However, the time of dispersal is highly variable and it is thus not unlikely that fully adult sons may at times be captured along with their fathers when multiple groups are sampled. Moreover, at least in some species such as F. micklemi, not all breeding males develop perioral stains. This complicates the identification of the breeding male based on perioral secretions alone and calls for a careful diagnosis of breeding status that takes other traits into account. Besides perioral stain, various studies report that male breeding status can be assessed by palpation of the testes, which are assumed to be larger in breeding males (e.g. 8 ,). However, available data on whether relative testes mass is greater in reproductive compared to non-reproductive Fukomys males are conflicting 31,58 .
In case of doubt, a promising indicator might instead be the width of the upper incisors. Observations on captive mole-rat families suggest that incisor width is a good relative age marker, particularly in males ( 59 ; pers. obs.). For Zaisan mole voles (Ellobius tancrei), a subterranean murid species with bathyergid-like extrabuccal incisors, it has already been demonstrated that incisor width can act as an age marker well into adulthood, before www.nature.com/scientificreports/ values for different age cohorts will eventually converge 60 . At what age incisor width becomes uninformative to differentiate breeding males from younger non-breeders in Fukomys remains to be determined.

Olfactory preference tests.
As indicated by non-significant differences in total sniffing time across experimental conditions, the mole-rats' interest in foreign conspecific odors was comparable across the three test situations. However, dependent on the offered scents, subjects preferentially allocated the time spent sniffing to one of the two options. Male perioral secretions were found to be preferred by conspecifics over scent samples taken from the same individual's dorsum and to exhibit a comparable attractiveness to anogenital odor. Anogenital smear evidently conveys rich social information in Fukomys 13,16 and exceeds other body scents in its quality to effectively signal individual identity 43 . The equivalence of perioral and anogenital swabs in the preference assay thus suggests an important communicative role for male perioral secretions in Fukomys. For female perioral swabs, we did not find a significant preference over samples taken from the back pelage. This could indicate a difference in the perceived informative value of female perioral odor or might simply be the result of less secretions being produced by females, resulting in a fainter scent. A possible bias in our behavioral assays might be that we have not considered the reproductive status of neither the donor nor the test subjects. Our subsequently generated results on perioral metabolomics suggest that mole-rats could be able to differentiate these status groups based on sebum volatiles and therefore might adopt their behavior accordingly (see below). However, as we quantified the relative informative value of odor sample pairs derived from the same respective individuals, we would not expect this issue to be a notable confounding factor here (see also 61 ). Nevertheless, the reproductive status of scent-sampled individuals should definitely be considered in future studies on these animals.
In any case, these results from the assays further indicate an asymmetric signaling function for perioral odors, with males investing more in the quantity of secretions to convey socially relevant signals to conspecific receivers than females do. The observation that females spent notably longer examining male secretions than the opposite sex did, might suggest that male perioral sebum has a particular role in intersexual communication.

Metabolomic profiling. Our analyses of volatile compounds detected via GCxGC-MS demonstrate notable
individual variation and striking sexual dimorphism in the volatile chemical composition of perioral secretions in Fukomys. Furthermore, they provide preliminary evidence for reproductive status-dependent signaling in both sexes. However, greater sample sizes are needed to robustly confirm this pattern and to confidently identify compounds that differentiate between intrasexual reproductive status groups, particularly male ones. Although sexual signatures were highly significant, it should be pointed out that the majority of compounds, in particular the ones occurring in the highest concentrations, are found among both sexes and can thus be expected to signal species identity. Highly sex-specific compounds represented only a fraction of the diversity and quantity of detected volatiles. Therefore, the proportional mixture of several compounds likely conveys a sexual signal in Fukomys perioral secretions, a pattern also known from the scent glands of various other rodents 62 .
Our GCxGC-MS approach used a method which compares mass spectra and retention indices of the detected compounds from perioral secretions to those in an existing library to identify volatiles. Compound names are thus just the most likely estimates. However, some of these compounds have been intensively studied and represent 'good matches' even without using external standards. Many of the metabolites that we detected were previously characterized in various other organisms, including bacteria (see below). This suggests that important fractions of the recovered compound diversity are not of endogenous origin but are produced by microbes colonizing the perioral sebum. This aligns well with the finding that sebaceous secretions in other mammals also house rich microbiota 63 . The notably small overlap in compounds between sampled food items and sebum demonstrates that food-derived odorants did not significantly bias our analyses. Nevertheless, besides compounds of endogenous and microbial origin, odorants presented in the perioral area might also derive from yet other sources. For instance, the sebum might act as a hydrophobic sponge that adsorbs additional odoriferous compounds from urine or feces that are transferred to the perioral region during (auto)coprophagy or anogenital autogrooming. Future studies should aim to clarify the exact origins of odorous compounds from the perioral sebum and what information mole-rats can deduce from them.
The most abundant volatile in both sexes was para-Cresol which conveys a typical 'pig smell' and is also secreted by male elephants during musth 64 . We also detected two variants of Bisabolenes (Bisabolene < (E)-gamma-> and beta) which were highly abundant. These sesquiterpenes are known to be produced by many plants, as well as by microbes and fungi, and to act as pheromones in insects 65 . Caryophyllene (3rd most abundant compound) is a natural sesquiterpene, present for example in cannabis and hops, which has a high affinity to the CB2 receptor in mice, with strong anti-inflammatory effects 66 . Similarly, 2-Heptanone is a ketone that stimulates alarm reactions in insects, while it evokes anxiety reactions in mice and rats, even without involvement of the vomeronasal organ 67 . This compound was also found to be abundant in both sexes. Significantly sex-biased metabolites also showed great structural diversity. For instance, tetramethyl-Pyrazine is a bacterial metabolite and significantly female-biased in our data. Likewise, 1-Phenyl-2-butanone is significantly female-biased. It was previously detected in defensive secretions of various invertebrates including millipedes 68 . Sulcatone is also female-biased and represents a ubiquitous eukaryote metabolite. 5-Hexadecanolide is significantly male-biased in our data and is known to otherwise act as a pheromone in the queens of the Oriental hornet (Vespa orientalis) 69 . 1-Ethenylaziridine is male-biased as well and has previously been detected in various bacterial species 70 .
An influence of reproductive and/or social status on volatile compounds of sebaceous gland secretions, as we observed in giant mole-rats, has also been demonstrated in a number of other social mammals, including rodents 71,72 , primates 73 , and carnivorans 63 . Respective odor profiles are hypothesized to provide an honest signal of rank and physiological condition to potential competitors 72 www.nature.com/scientificreports/ aggression. Sex and status-dependent signals from the perioral glands might serve this role in Fukomys families as well, facilitating the identification and social evaluation of both group members and foreign individuals that might enter an established family (compare 61 ). Whether or not differences in the odor profiles between reproductive status groups are adaptive, they might proximately be determined by the hormonal status of a respective individual. It has been shown that sex hormone levels as well as pregnancy and lactation can affect commensal microbial communities and do regulate body odor via this path 71,73,74 . To which extent these factors might explain differentiation in Fukomys perioral volatile diversity remains to be determined. Given that perioral nuzzling is typically initiating copulation 29 , it might be tempting to speculate that certain volatiles act as sex pheromones in mole-rats. Indeed, some of the compounds we detected have been suggested to represent sexual pheromones in mice, for instance 2-Heptanone 75 . However, contrary to expectation, these molecules show no sexually dimorphic expression pattern in giant mole-rats. Apart from that, it should be noted that the vomeronasal system, which typically responds to pheromones, appears to be of little relevance in Fukomys and other African mole-rats. The bathyergid vomeronasal organ is growth-deficient 76,77 and its vomeronasal receptor repertoire is small. The latter is typical for subterranean rodents in general 78 . Nevertheless, there are still structural indications for the bathyergid vomeronasal organ being functional 76 and even if that is not the case, pheromones may also be perceived by the primary olfactory epithelium 79 . Hence, as of right now, a pheromone function for perioral compounds cannot be ruled out.

Synopsis and conclusion
We have shown that Northern common mole-rats (as well as naked mole-rats) of either sex possess complex and structurally derived sebaceous glands situated in their mouth corners. These perioral glands display a sexspecific secretion activity, which is higher in males and not affected by reproductive status. Conspecifics of both sexes, but particularly females, are notably responsive to male perioral secretion, which suggests them serving in social communication. Finally, metabolomic profiling revealed that the composition of volatile compounds from perioral sebum is sexually dimorphic and also allows to distinguish breeders and non-breeders.
These results suggest that perioral secretions convey important social information in both sexes but also that the investment into perioral signaling is asymmetrical. It appears possible that while female secretions might be exclusively perceived during close contact (i.e. facial nuzzling) the greater quantities of male secretion enable a more potent olfactory signal that might play a role in attracting and courting mates. It is tempting to speculate that scents derived from the perioral glands are passively deposited onto the soil during tooth digging. The characteristic waxy consistency of the secretion might aid in prolonging the longevity of the scent signal (compare 80 ). This way, for instance, male tenants could effectively signal their presence in a burrow system and perhaps also their physical condition to same-sex intruders that might challenge their position 5,8 .
In any case, the perioral secretions of Fukomys can be added to a long list of sexually dimorphic traits in these monogamous, cooperatively-breeding rodents, all pointing to a notable role of male intrasexual competition within their social system 11 .
Several novel questions on olfactory signaling in common mole-rats and other bathyergids arise from this research. For instance, it remains to be determined how perioral and anal gland derived scents complement each other in the mediation of social behaviors and whether traces of perioral sebum can indeed act as lasting scent marks in burrow systems. The intriguing intrasexual variability in the expression of perioral stains, together with its physiological causes and behavioral implications, also calls for further clarification. It is to hope that such work on the conspicuous perioral secretions of social mole-rats will not only clarify how these animals effectively communicate underground, but also stimulate research on the understudied mouth corner glands of other rodent species.

Data availability
The paper and its accompanying supplementary information contain all data discussed in the study. www.nature.com/scientificreports/ Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.